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I. INTRODUCTION 

Quantum tunneling, where a particle has finite prob- 
ability to penetrate a classically forbidden region is an 
important feature of wave mechanics. Invention of the 
tunnel diode 0, the scanning tunneling microscope 
etc. make it useful from a technological point of view. In 
1932 MacColl pointed out that tunneling is not only 
characterized by a tunneling probability but also by a 
time the tunneling particle takes to traverse the barrier. 
There is considerable interest on the question of time 
spent by a particle in a given region of space 0, 0, . 
The recent development of nanotechnology brought new 
urgency to study the tunneling time as it is directly re- 
lated to the maximum attainable speed of nanoscale elec- 
tronic devices. In a number of numerical 0, experimen- 
tal [H, and analytical study of quantum tunneling pro- 
cesses, various definitions of tunneling times have been 
investigated. These different time scales are based on 
various different operational definitions and physical in- 
terpretations. Till date there is no clear consensus about 
the existence of a simple expression for this time as there 
is no hermitian operator associated with it 0. Among 
the various time scales, 'dwell time' 01 which gives the 
duration of a particle's stay in the barrier region regard- 
less of how it escapes can be calculated as the total prob- 
ability of the particle inside the barrier divided by the 
incidentprobability current. Biittiker and Landauer pro- 
posed [Tl| that one should study 'tunneling time' using 
the transmission coefficient through a static barrier of in- 
terest, supplemented by a small oscillatory perturbation. 
A large number of researchers interpret the phase delay 
time j^, IT2I IT^ I as the temporal delay of a transmitted 
wave packet. This time is usually taken as the differ- 
ence between the time at which the peak of the transmit- 
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ted packet leaves the barrier and the time at which the 
peak of the incident Gaussian wave packet arrives at the 
barrier. Within the stationary phase approximation the 
phase time can be calculated from the energy derivative 
of the 'phase shift' in the transmitted or reflected am- 
plitudes. Buttiker-Landauer raised objection that 
the peak is not a reliable characteristic of packets dis- 
torted during the tunneling process. In contrast to 'dwell 
time' which can be defined locally, the 'phase time' is es- 
sentially asymptotic in character The 'phase time' 
statistics is intimately connected with dynamic admit- 
tance of micro-structures 0. This 'phase time' is also 
directly related to the density of states 0] . The univer- 
sality of 'phase time' distributions in random and chaotic 
systems has already been established earlier [l7|. In the 
case of 'not too opaque' barriers, the tunneling time eval- 
uated either as a sample 'phase time' or calculated 
through the analysis of the wave packet behaviour 0] 
becomes independent of the barrier width. This phe- 
nomenon is termed as the Hartman effect ITcl Il9|. 
This implies that for sufficiently large barriers the effec- 
tive velocity of the particle can become arbitrarily large, 
even larger than the light speed in the vacuum (super- 
luminal effect). Though this interpretation is a little far 
fetched for non-relativistic Schrdinger equation as veloc- 
ity of light plays no role in it, this effect has been estab- 
lished even in relativistic quantum mechanics. 

Though experiments with electrons for verifying this 
prediction are yet to be done, the formal identity be- 
tween the Schrodinger equation and the Hclmholtz equa- 
tion for electromagnetic wave enables one to correlate the 
results for electromagnetic and microwaves to that for 
electrons. Photonic experiments show that electromag- 
netic pulses travel with group velocities in excess of the 
speed of light in vacuum as they tunnel through a con- 
striction in a waveguide |2f| . Experiments with photonic 
band-gap structures clearly demonstrate that 'tunneling 
photons' indeed travel with superluminal group veloci- 
ties Q. Their measured tunneling time is practically 
obtained by comparing the two peaks of the incident and 



transmitted wave packets. Thus all these experiments di- 
rectly or indirectly confirmed the occurrence of Hartman 
effect without violating 'Einstein causality' i.e., the sig- 
nal velocity or the information transfer velocity is always 
bounded by the velocity of light. This tunneling time 
could be interpreted as the 'time of passage' of the peak. 
Since the velocity of the 'peak' may exceed arbitrarily 
large numbers, this 'fast tunneling' has been frequently 
related to 'superluminal propagation' |2£j. The 'Hart- 
man effect' has been extensively studied both for nonrel- 
ativistic (Schrodinger equation) and relativistic (Dirac 
equation) [il IH l20j| cases. Recently Winful snowed 
that the saturation of phase time is a direct consequence 
of saturation of integrated probability density under the 
barrier. Due to this saturation addition of a new particle 
in the incident side leads to an almost immediate release 
of another particle from the other side and these two par- 
ticles are causally unrelated; i.e. even the superluminal 
tunneling does not violate causality. Equivalently, for 
electromagnetic waves the origin of the Hartman effect 
has been traced to stored energy. Since the stored energy 
in the evanescent field decreases exponentially within the 
barrier after a certain decay distances it becomes inde- 
pendent of the width of the barrier. The Hartman effect 
has been found in one dimensional barrier tunneling [T^ | 
as well as for cases beyond one dimension as in tunneling 
through mesoscopic rings in presence of Aharonov-Bohm 
flux |22( . In the current note we extend the study of phase 
times for branched networks of quantum wires. Such net- 
works can readily be realized in optical wave propagation 
experiments. 

In Section-HTl we give the description of typical quan- 
tum network systems under consideration. The theoreti- 
cal framework used to analyze these systems is provided 
m Section-lrm All the main results showing Hartman ef- 
fect and the effects of quantum non-locality are discussed 
m 

Section-El Finally we summarize the results and draw 
conclusions in Section-Tvl 

II. DESCRIPTION OF THE SYSTEM 

As a model system, we choose a network of thin wires. 
The width of these wires are so narrow that only the 
motion along the length of the wires is of interest (a single 
channel case) . The motion in the perpendicular direction 
is frozen in the lowest transverse subband. In a three-port 
Y-branch circuit (Fig. ^) two side branches of quantum 
wire 5*i and 5*2 are connected to a 'base' arm So at the 
junction J. In general one can have iV(> 2) such side 
branches connected to a 'base' wire. 

We study the scattering problem across a network ge- 
ometry as presented schematically in Fig. ^ 

Such geometries are important from the point of view 
of basic science due to their properties of tunneling and 
interference [23l |2^ | as well as in applications such as 
wiring in nano-structures. In particular, the Y-junction 
carbon nanotubes are in extensive studies and they show 
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FIG. 1: Schematic diagram of a Y-junction or three-way split- 
ter. 



various interesting pro per ties like asymmetric current 
voltage characteristics 25]. In our system of interest 
there are finite quantum mechanical potential barriers of 
strength V n and width w n in the n-th side branch. The 
number of side branches can vary from n = 2, 3, • • ■ N. 
We focus on a situation wherein the incident electrons 
have an energy E less than V n for all n. The impinging 
electrons in this subbarrier regime travels as an evanes- 
cent mode/ wave and the transmission involve contribu- 
tions from quantum tunneling and multiple reflections 
between each pair of barriers and the junction point. 
Here we are interested in a single channel case where 
the Fermi energy lie in the lowest subband. To excite the 
evanescent modes in the side branches one has to produce 
constrictions by making the width of the regions of wires 
containing barriers much thinner than that of other parts 
of the wires. The electrons occupying the lowest subband 
in the connecting wire on entering the constrictions ex- 
perience a potential barrier (due to higher quantum zero 
point energy) and propagate as an evanescent mode [27| ■ 
In this work an analysis of the phase time or the group 
delay time in such a system is carried out. 



III. THEORETICAL TREATMENT 

We approach this scattering problem using the quan- 
tum wave guide theory [2|| [23] . In the stationary case 
the incoming particles are represented by a plane wave 
e ikx Q f un ^ amplitude. The effective mass of the prop- 
agating particle is m and the energy is E = h 2 k 2 /2m 
where k is the wave vector corresponding to the free 
particle. The wave functions, which are solutions of the 
Schrodinger equation, in different regions of the system 
considered (Fig. can be written as, 
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^n(zo) = e ik *° + Re- ikx ° (mS ), 

*P(n)A x n) = A n e lkXn + B n e' lkXn (region I in S n ) , 

ip {n)u {x n ) = C n e- K ^ x -- lK) + £> n e«» (region II in 5„) , 

TpinhtAxn) = t n e ik ^- lb «- w ^ (region III in S n ) , 



with K n — y 2m(V n — E) /h 2 being the imaginary wave 
vector in presence of rectangular barrier of strength V n . 
'4>(n) I , ^(n)// an d ^(njm denote wave functions in three 
regions /, II and III, respectively, on n-th side branch. 
xq is the spatial coordinate for the 'base' wire, whereas 
x n is spatial coordinate for the n-th arm. All these coor- 
dinates are measured from the junction J. In rt-th side 
branch, the barrier starts at a distance lb n from the junc- 
tion J. 

We use Griffith's boundary conditions [2^| 



4>in{J) = ^(n=l)i( J ) = ~4>{n=2) 1 {J) 

and 

dlpin{x ) 



dx 



— — 7{ 

j ox n 



■>j\n=N) I (J), 
(1) 



(2) 



at the junction J. All the derivatives are taken either 
outward or inward from the junction |26j . In each side 
branch, at the starting and end points of the barrier, the 
boundary conditions can be written as 



■0(n)„ {IK + W n ) = t/}( n ) III (lb 



dXn 



dXr 



(ib n ) dx 



(lb n ) 



dx r , 



{lb n +w rl ) 



(3) 
(4) 

(5) 
(6) 



From the above mentioned boundary conditions one 
can obtain the complex transmission amplitudes t n on 
each of the side branches. 



IV. RESULTS AND DISCUSSIONS 

Once t n is known, the 'phase time' (phase time for 
transmission) can be calculated from the energy deriva- 
tive of the phase of the transmission amplitude 0, 0] as 



In what follows, let us set U = 1 and 2m — 1. We now 
proceed to analyze the behavior of r„ as a function of 
various physical parameters for different network topolo- 
gies. We measure time at the far end of each barrier in 
the branched arms containing barriers and in the case 
of arms in absence of any barrier we measure the phase 
time at the junction points. We express all the physi- 
cal quantities in dimcnsionless units i.e. all the barrier 
strengths V n in units of incident energy E (V n = V n /E), 
all the barrier widths w n in units of inverse wave vector 
k~ 1 (w n = kw n ), where k = \f~E~ and all the extrapo- 
lated phase time t„ in units of inverse of incident energy 
E {r n = Er n ). 

First we take up a system similar to the Y-junction 
shown in Fig. ^ in presence of a barrier V\ of width w\ 
in arm S\ but in absence of any barrier in arm S*2. For 
a tunneling particle having energy E < V\ we find out 
the phase time T\ in arm S\ as well as in arm Si as a 
function of barrier width w\ (Fig. From Fig. [2fa) it 
is clear that n evolves with w\ and eventually saturates 
to t s i for large w\ to show the Hartman effect. Fig. Etb) 



(a) 4 -5 

4 

3.5 

H 3 
2.5 
2 
1.5 



(b) 

2.5 

2 

H 1.5 I 
1 

0.5 I 
0- 



5 10 15 20 25 30 



5 10 15 20 25 30 



dArg[t v 
dE 



(7) 



where, v = hk/m is the velocity of the free particle. 
The concept of 'phase delay time' was first introduced 
by Wigner [T^ | to estimate how long a quantum mechan- 
ical wave packet is delayed by the scattering obstacle. 



FIG. 2: For a 3-way splitter with a barrier in Si arm, the 
'phase times' n and T2 are plotted as a function of barrier 
width 'wi' in (a) and (b) respectively. The solid, dotted, 
dashed, dot-dashed and the dashed-double dotted curves are 
for Vi = 5, 4, 3, 2, and 1.05 respectively. Other system param- 
eters are E = 1, Ibi = 3. 
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FIG. 3: Here for a 3- way splitter with one barrier in each 
branched arm Si and S2, the 'phase times' ti is plotted as a 
function of barrier width 'wi' keeping W2(= 1) and Vi(= 5) 
fixed and for different values of parameter V2. The small 
dashed, dot-dashed, solid, long dashed and dotted curves are 
for V% = 2.5, 3.5, 5.0, 10.0 and 12.5 respectively. Other system 
parameters are E = 1, Ibi = Z&2 = 3. In the inset r s i is plotted 
as a function of V2 for the same system parameters. 
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FIG. 4: Here for a 3- way splitter with one barrier in each side 
branch Si and 5*2, the 'phase times' n (solid curve) and \ti\ 2 
(dashed) are plotted as a function of 'wi' for a very small 
W2(— 0.5). Other system parameters are E — 1, lb\ — W2 = 
2.5, V2 = 5 and Vi = 15. In the inset, the solid and dashed 
curves represent T2 and \t2\ 2 respectively as a function of wi. 
For better visibility we have plotted phase times scaled down 
by a factor of 20. 



shows the phase time T2 in arm S2 which does not contain 
any barrier. This also evolves and saturates with w\, the 
length of the barrier in the other arm S\. This delay is 
due to the contribution from paths which undergo multi- 
ple reflection in the first branch before entering the sec- 
ond branch via junction point J. In absence of a barrier 
in the n-th arm the phase time T n measured close to the 
junction J should go to zero i.e. r„ — ► in the absence of 
multiple scatterings in the first arm. Note that t s \ and 
t S 2 change with energies of the incident particle (Fig.[3J). 
From Fig|21it can be easily seen that t S 2 is always smaller 
than T s i for any particular Vi i. e. the saturation time in 
the arm having no barrier is smaller. The phase time in 
both the arms show non-monotonic behavior as a func- 
tion of V\. As we decrease the strength of the barrier 
V\ the value of t\ fa) decreases in the whole range of 
widths of the barrier and also the saturated value of 77,1 
(t S 2) decreases until V\ reaches 1.6 and with further de- 
crease in V\ the values of t\ fa) as well as t s i (t S 2) starts 
increasing. 

As the second case we take up another Y-j unction 
which contain potential barriers in both its side branches 
as shown in Fig.^ We fix the values of Vi(= 5) and vary 
wi for each values of V2 to study the wi-dependence of 
Ti (Fig|3J. From Fig. [3] we see that t% decreases with 
increase in w± to saturate to a value r s i at each value of 
V2 thereby showing 'Hartman effect' for arm 'Si'. But 
now, we can tune the saturation phase time at arm Si 
non-locally by tuning strength of the barrier potential 
V2 sitting on another arm S2! Thus 'quantum nonlocal- 
ity' enables us to control the 'super arrival' time in one 
of the arms (Si) by changing a parameter (V2) on an- 



other, spatially separated from it. In the inset of Fig. |3| 
we plot t s i as a function of V2. It clearly shows that when 
the barrier strengths Vi and V2 are very close the 'phase 
time' reaches its minimum value. In all other cases i.e. 
whenever Vi 7^ V2, the value of r s i is larger. 

We will show now another interesting result related 
to the Hartman effect. For this we keep V^(= 5) un- 
altered and reduce u>2- For very small W2{= 0.5) we see 
from Fig.0]that n is negative for almost the whole range 
of uii-values showing 'time- advancement' and eventually 
after a sharp decrease saturates to a negative value of 
T a i = —4.514 implying 'Hartman effect' with advanced 
time. It might be noted that, in principle, the 'time- 
advancement' (Fig. ^ can be measured experimentally 
as |ti| 2 has a non-zero finite value for a small range of 
i«i at lower wi regime where n is negative. In the inset 
we plot the corresponding T2 and | £2 1 2 as function of w\. 
Again the values of t~2 remains different from the one 
dimensional tunneling through a barrier of strength V2 
and width W2 in the whole range of wi implying 'quan- 
tum nonlocality'. In the cases discussed so far T2 vary 
more sharply in small wi regime. Further the inset in 
Fig.0|shows a dip in T2 at parameter regimes where \t2\ 2 
has a minimum. For a wave packet with large spread 
in real space it is possible that the leading edge of the 
wave packet reaches the barrier much earlier than the 
peak of the packet. This leading edge in turn can tun- 
nel through to produce a peak in the other end of the 
barrier much before the incident wave packet reaches the 
barrier region, sometimes referred to as pulse reshaping 
effect. This, in general, causes 'time advancement 
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FIG. 5: Here iV(> 2) number of side branches each with a 
barrier of strength V(= 5) and width w(= 100) are connected 
with the incident arm So- In each arm barriers start at the 
same point lb{= 1). Thus all the side branches being identical 
the 'phase time' for transmission through these arms S n , n = 
2, 3, ■ ■ ■ N saturate to the same value t s . The saturated 'phase 
time' t s is plotted as a function of the total number of side 
branches N in the system. The dotted and solid curves are 
for V = 5.0 and 1.25 respectively. The incident energy is kept 
at E = 1. 



This negative delay does not violate causality, however, 
the time is bounded from the below. In the presence 
of square wells in one dimensional systems negative time 
delays have been observed. This effect is termed as 'ultra 
Hartman effect' [ see for details |30J ]. 

Finally consider a similar system as that shown in 
Fig. ^ but in presence of N(> 2) identical side branches 
and study phase time as a function of increasing N. 
All the side branches being identical the 'phase times' 
for transmission through each of these arms S n ,n = 
2,3,- • ■ N saturate to the same value t s for very large 
w n . In Fig|S]we plot the saturation value r s as a func- 
tion of the total number of side branches N present in 
the system. From the figure we see that for V = 5, t s 
first increases with N to a maximum value of 3.776 at 
N = 9 and thereafter keeps on decreasing with the in- 
crease of N. As we start reducing the strength of the 
barriers from 5 we see that for V = 1.49 the increasing 
nature of r s in small N range vanishes. In general, at 
larger N, the decreasing nature of t s with N persists, 
e.g., note the solid curve in Fig[S] plotted at V — 1.25, 
but the initial increase in r s is not a generic feature. For 
larger N transmission amplitude in each side branch re- 



duces with increase in N and hence the corresponding 
peaks of wave packets reach the far end at earlier times 
thereby reducing r s . 

V. CONCLUSIONS 

We have studied Hartman effect and non-locality in 
quantum network consisting of a main one dimensional 
arm having 7V(> 2) side branches. These side branches 
may or may not have barriers. In presence of barrier 
the 'phase time' for transmission through a side branch 
shows the 'Hartman effect'. In general, as the num- 
ber of side branches N increases, the saturated 'phase 
time' decreases. Due to quantum nonlocality the 'phase 
time' and it's saturated value at any side branch feels 
the presence of barriers in other branches. Thus one 
can tune the saturation value of 'phase time' and conse- 
quently the superluminal speed in one branch by chang- 
ing barrier strength or width in any other branch, spa- 
tially separated from the former. Moreover Hartman ef- 
fect with negative saturation times (time advancement) 
has been observed for some cases. In conclusion gen- 
eralization of Hartman effect in branched networks ex- 
hibits several counter-intuitive results due to quantum 
non-locality. System parameters such as number of side 
branches N and barrier widths w n , strengths V n , distance 
lb n (from J) and incident energy E etc. play very sen- 
sitive roles in determining delay times. The delay times 
are also sensitive to the junction S-matrix elements used 
for a given problem. In our present problem junction S'- 
matrix is determined uniquely by the wave guide trans- 
port methods. Depending on lb n there may be one or 
several bound states located between the barriers in dif- 
ferent branched arms and as a consequence saturated de- 
lay time can be varied from the negative (ultra Hartman 
effect) to positive and vice-versa. We have verified this 
by looking at the transmission coefficient in the second 
arm S2 which exhibits clear resonances as a function of 
Z62 ■ Around these resonances, saturated delay time in 
the first arm r s i changes the sign [the details of which 
will be published elsewhere H2]- Moreover the reported 
effects are amenable to experimental verifications in the 
electromagnetic wave-guide networks. 
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